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ABSTRACT 

IC 10 X-l is a variable X-ray source i n the Local Group sta rburst galaxy IC 10 whose optical 
counterpart is a Wolf-Rayet (WR) star. iPrestwich et alJ (|2007f ) recently proposed that it contains 
the most massive known stellar-mass black hole (23-34 M©), but their conclusion was based on radial 
velocities derived from only a few optical spectra, the most important of which was seriously affected 
by a CCD defect. Here we present new spectra of the WR star, spanning one month, obtained with the 
Keck-1 10 m telescope. The spectra show a periodic shift in the He II A4686 emission line as compared 
with IC 10 nebular lines such as [O III] A5007. From this, we c alculate a period of 3 4.93 ± 0.04 hr 
(consistent with the X-ray period of 34.40 ± 0.83 hr reported bv IPrestwich et aT1l2007l ) and a radial- 
velocity semi - amplitude of 370 ±20 km s -1 . The resulting mass function is 7.64 ±1.26 M©, consistent 
with that of IPrestwich et ail (|2007f ) (7.8 Mq). This, combined with the previously estimated (from 
spectra) mass of 35 M© for the WR star, yields a minimum primary mass of 32.7 ± 2.6 M Q . Even if 
the WR star has a mass of only 17 M©, the minimum primary mass is 23.1±2.1 M©. Thus, IC 10 X-l 
is indeed a WR/black-hole binary containing the most massive known stellar-mass black hole. 
Subject headings: galaxies — starburst, stars — Wolf-Rayet, X-rays — binaries, black hole physics 



1. INTRODUCTION 

IC 10 X-l is a bright, variable X-ray source in the Lo- 
cal Group metal-poor starburst ga laxy IC 10 w i th an 
X-ray luminosity of 1 38 erg s" 1 (IBrandt et all 119971 : 
iBauer fc Brandtj l2004h . lLozinskava fc Moiseevl (|2007l ) 
suggested that IC 10 X-l is the compact remnant of 
a hypernova, based on the nature of the synchrotron 
superbubble in IC 10. There are four possible optical 
counterparts to the X-ray source, the most likely be- 
ing the luminous Wol f-Rayet (WR) star [MAC92] 17A 
(|Crowther et al.ll2003l ). Previous spectroscopic observa- 
tions of [MAC92] 17A s h ow pr ominent He II line emis- 
sion; IClark fc Crowtherl (|2004f ) classified it as a WNE 
star. 

IPrestwich et alJ (|2007l ) recently proposed that IC 10 X- 
1 and [MAC92] 17A are a WR/black-hole (BH) binary 
containing the most massive known stellar-mass black 
hole (23-34 M©). Unfortunately, their conclusion was 
based on radial velocity measurements from only a few 
optical spectra, and the observation showing an apparent 
spectral shift was seriously affected by a CCD defect, 
casting some doubt on the reality of the shift. They also 
assumed that the observed X-ray period of IC 10 X-l 
(34.40 ± 0.83 hr) is equal to the orbital period of the 
binary system. 

Here we present radial velocities from 10 new optical 
spectra of the WR star sp anning one month; our results 
confirm the conclusions of IPrestwich et all (|2007f ). In 
we describe the observations and data reduction, in [J3]we 
discuss our analysis of the spectra, and in f}4] we present 
ou r results. A preliminary report on this work is given 
bv lSilverman fc Filippenkol (|2008[ ). 

2. OBSERVATIONS AND DATA REDUCTION 

[MAC92] 17A was observed with the Low R esolution 
Imaging Spectrometer (LRIS; lOke et aT1ll995h mounted 
on the Keck-I 10 m telescope in 2007 during the nights 
of Nov. 11-13 (UT dates are used throughout this pa- 



per), Nov. 16, and Dec. 12; a journal of observations is 
given in Table [TJ LRIS is now equipped with an atmo- 
spheri c dispersion com pensator; thus, differential light 
losses (|Filippenkdll982l ) were not a problem, even at high 
airmasses. 

Fortuitously, th e long slit included R SMV 2, another 
WR star in IC 10 (jCrowther et al.ll2003l ). Unfortunately, 
the light from this second object fell between the two 
chips on the blue side of LRIS in the observations on 
Nov. 11.492 and Nov. 13.320. 

All data were taken using a 5600 A dichroic and a 
600/4000 grism on the blue side of the spectrograph. All 
but one observation employed a long slit of width 0'.'7; 
the observation on 2007 Nov. 16.242 used a long slit of 
width l'.'O. The slit was always oriented at a position 
angle (PA) of 162°, through another star readily visible 
on the acquisition and guider image, in order to facilitate 
object acquisition and to minimize contamination from 
adjacent stars. The spectra span a wavelength range 
of ~ 3800-5550 A, and all observations had a FWHM 
spectral resolution of ~3.5 A (^220 km s _1 ). 

The data were r educed using standard techniques (e.g., 
IFolev et al.l f2003h . Routine CCD processing and spec- 
trum extraction for the data were completed with IRAF 1 . 
T he data were extracted with the optimal algorithm 
of I Home] (|1986f ). We obtained the overall wavelength 
scale from low-order polynomial fits to calibration-lamp 
spectra. Small wavelength shifts were applied to in- 
dividual spectra after cross-correlating night-sky lines 
with a template sky spectrum. Using our own IDL 
routines, we fit spectrophotometric standard-star spec- 
tra to flux-calibrate our data and remove telluric lines 
(|Wade fc Horndfl988t iMatheson et al.ll2000f) . 

1 IRAF, the Image Reduction and Analysis Facility, is dis- 
tributed by the National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc. (AURA) under cooperative agreement with the Na- 
tional Science Foundation (NSF). 
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TABLE 1 

Journal of Observations and Radial Velocities 



HJD a 


2007 UT Date b 


Exp. (s) 


Airmass c 


Seeing 


Phase 


Me 11 v (km s j 


TT TT DU/IHf / A \ 

He 11 r WrlM (A) 


2454415.865 


Nov. 


11.361 


1400 


1.33 


0.8 


0.9801 


384 ± 29 


14.3 ±0.9 


2454415.882 


Nov. 


11.378 


1400 


1.36 


0.8 


0.9918 


362 ± 24 


13.3 ±0.8 


2454415.996 


Nov. 


11.492 


1800 


1.95 


1.5 


0.0701 


325 ± 40 


14.1 ± 1.5 


2454416.761 


Nov. 


12.257 


1800 


1.35 


0.7 


0.5957 


-312 ± 44 


18.1 ± 1.7 


2454416.932 


Nov. 


12.428 


1800 


1.53 


0.8 


0.7132 


-85 ±43 


16.3 ±1.7 


2454417.824 


Nov. 


13.320 


1800 


1.30 


1.5 


0.3260 


-182 ±26 


15.5 ±0.9 


2454420.708 


Nov. 


16.205 


1500 


1.46 


1.0 


0.3081 


-77 ± 113 


19.7 ±4.6 


2454420.726 


Nov. 


16.223 


1500 


1.40 


1.1 


0.3205 


-148 ±52 


13.6 ±2.0 


2454420.745 


Nov. 


16.2428 


1500 


1.36 


1.0 


0.3335 


-133 ±51 


14.0 ±1.9 


2454446.761 


Dec. 


12.258 


1200 


1.31 


1.4 


0.2074 


91 ±51 


24.6 ±2.1 



Airmass at midpoint of 
34.93 ± 0.04 hr and T 




a Heliocentric Julian Date at midpoint of exposure. UT Date at midpoint of exposure.' 
exposure. d Approximate full width at half-maximum intensity (FWHM), arcsec. c Using P = 

2007 Nov. 9.935 ± 0.014 (the UT date of maximum radial velocity).' Radial velocity calculated relative to the [O III] 
A5007 line. 8 This is the only observation with a long slit of width l'.'O; all others used a slit of width 0'.'7. 

formed in the inner part of the WR wind (close to the 
star), justifying i ts use when determini ng radial velocities 
of WR binaries (|Prestwich et al.ll2007f ). 

For each observation we fit a Gaussian profile to the 
nebular [O III] and H/3 lines, as well as to the He II line of 
[MAC92] 17A. To remove possible errors introduced by 
our wavelength and flux calibrations, the profiles were 
fit to the raw, one-dimensional (1-D) extracted spectra 
(i.e., neither wavelength nor flux calibrated). Owing 
to the small number of emission lines available in the 
calibration-lamp spectra, the wavelength solution was 
found to vary quite a bit from spectrum to spectrum, 
when in fact nothing was actually changing except per- 
haps the zero point. Fitting profiles to this raw 1-D data 
is reasonable for our purposes because the radial velocity 
of the WR star only depends on the relative change in 
the spacing between the He II line and the (stationary) 
nebular lines. 

Figure shows three spectra, centered on the He II 
A4686 spectral feature, in velocity space (as calculated 
relative to the [O III] A5007 line). The signal-to-noise 
ratio (S/N) of the observations is typical of our entire 
dataset. It is clear that the centroid of the He II line is 
shifted between the three observations. 

The FWHM of the He II line appears to change with 
time in Figure [H As seen in the last column of Table [TJ 
the only observation that actually has a substantially dif- 
ferent FWHM is the one fro m 2007 Dec. 12.2 5 8 (the mid- 
dle spectrum in Figured]). iPrestwich et alj (|2007l ) note 
that most of their epochs were well-fit by a Gaussian with 
a 15-17 A FWHM, which matches our data quite well (we 
have a mean of ~16.3 A). They also mention that one of 
their observations had a FWHM of only 10 A, but that 
the narrowness of this observation is likely due to a "CCD 
defect that contaminates the red component." None of 
our derived FWHM values are within ler of 10 A, suggest- 
ing that their assessment is co rrect. However, t his is also 
the only spectrum from which IPrestwich et alj (|2007f) de- 
duced a radial velocity variation in [MAC92] 17A, casting 
some doubt on the claimed shift. 

From our Gaussian fits we calculated the separation in 
pixels between the centroids of the He II feature and each 
host-galaxy nebular line. These were then converted to 
a difference in wavelength by using the dispersion value 
we derived from low-order polynomial fits to calibration- 
lamp spectra (0.621 A pixel -1 ). Next, we performed 



4000 4500 5000 5500 

Rest Wavelength (A) 

Fig. 1. — Unphased combined spectra of RSMV 2 (top) and 
[MAC92] 17A (bottom), and their superposed H II regions. Both 
nebular and WR line identifications are discussed in the text. The 
spectrum of RSMV 2 has been shifted up by 0.45 units for clarity. 

The unphased combined spectra of RSMV 2, from 8 
epoches, and of [MAC92] 17A, from all 10 epochs, are 
shown i n Figured) The radi al velocity of IC 10, — 348 ±1 
km s _1 (|Huchra et al.lfl999l ). has been removed. 

In the spectra of both objects, there are strong nebular 
emission lines of Hy, H/3, and [O III] (and US and He in 
RSVM 2) from H II regions in the slit. One can also see 
the relatively broad He II A4686 emission and the weaker 
He II A5411 emission which comes from the WR stars 
themselves. In addition, the fainter N V 4603-20 A and 
N III 4634-42 A emission features from [MAC92] 17A 
and RWVM 2 (respectively) are seen. Broad emission 
from relatively low-excitation species such as H I and 
He I is not present in [MAC92] 17A, suggesting that the 
emission does not arise from an accretion disk. 

3. ANALYSIS 

The X-ray period found bv IPrestwich et aTl (|2007h . as 
well as our measured period (see are substantially 

larger than the period needed for Roche lobe overflow, 
given reasonable masses of [MAC92] 17A and its compact 
binary companion. (We, however, derive an expected pe- 
riod for Roche lobe overflow of 9.0-15.9 hr, which differs 
from the 2-3.5 hr period suggested by Prestwich et al. 
2007). Thus, the accretion must be the result of a wind 
from the WR star, as is often seen in high-mass X-ray 
binaries. Furthermore, the He II A4686 emission line is 
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Fig. 2. — Spectra of the He II emission line from [MAC92] 17A 
obtained on 2007 Nov. 12.257 (top), 2007 Dec. 12.258 (middle), 
and 2007 Nov. 11.361 (bottom). The data have been deredshiftcd 
and converted into velocity space. The vertical long-dashed line is 
where the He II feature of [MAC92] 17A would be centered if at 
rest. The short-dashed lines show our Gaussian fits to the data. 
The centroid of the He II line clearly shifts with time. 

a nonlinear least-squares fit to the differences in wave- 
length. The resulting cosine yielded the period of the 
He II line and its systemic wavelength shift, which is 
equivalent to the zero point of the cosine. Finally, each 
spectrum's deviation from this systemic wavelength shift 
was converted into a radial velocity. Similarly, the semi- 
amplitude of the fit was converted from a wavelength to 
a velocity. 

4. RESULTS 
4.1. Radial- Velocity Curve 

Figure [3] shows the radial-velocity curve of 
[MAC92] 17A. Two complete periods, each with 10 
observations, are shown. The velocities in Figure [3] were 
calculated with respect to the [O III] A5007 line. We also 
determined radial velocities and an associated cosine 
fit with respect to the nebular H/3 line. These radial 
velocities differed from the ones derived using the [O III] 
line by much less than la at each epoch. In addition, 
the period of the curve derived using the H/3 line differed 
by only 0.01% from the curve shown in Figure [3j and 
the semi-amplitude differed by ~0.3%. These differences 
led to a mere 0.5% difference in the final BH mass. 

Since [O III] had a somewhat higher S/N than H/3 (as 
seen in Figuref!}, the H/3 Gaussian fits had slightly larger 
uncertainties in their centroids and thus led to larger 
overall uncertainties. Thus, we have chosen to present 
only the velocities derived relative to the [O III] line. 

The error bars shown represent ler, and were calculated 
by considering the uncertainty in the centroids of our 
Gaussian fits to both the He II line as well as the [O III] 
nebular line. The point in Figure [3] with the relatively 
large error bar is due to a low S/N spectrum. 
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Fig. 3. — Radial-velocity curve of [MAC92] 17A using velocities 
relative to the [O III] A5007 spectral line. Two cycles are shown for 
clarity. Formal velocity error bars arc la. Sec the text for values 
of the fit parameters. 
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Fig. 4. — Value of \ 2 versus trial period for [MAC92] 17A, using 
all 10 spectra discussed in the text. 

The four-parameter fit (zero point, semi-amplitude, 
period, and phase) yielded a semi-amplitude of K2 — 
370 ± 20 km s" 1 and a period of P = 34.93 ± 0.04 hr 
(consisten t with the X-ray perio d of 34.40 ± 0.83 hr re- 
ported bv lPrestwich et aLl l2007h However, the phasing 
of our radial- velocity curve with the X-ray light curve 
(|Prestwich et alj 120071 Fig. 1) is unknown because of 
the substantial uncertainty in the X-ray period. 

We examined the distribution of possible periods re- 
sulting from our observations in order to test whether 
our calculated period suffers from any ambiguity in the 
number of cycles between observations. Figure [4] shows 
X 2 versus trial period for our complete set of 10 obser- 
vations. It is clear that the two periods that best fit the 
data are ~34.9 hr ( X 2 « 2.43) and -35.7 hr ( X 2 « 2.50). 

We have chosen to adopt the 34.9 hr period since (a) 
it has the lowest value of x 2 (although the X 2 from the 
35.7 hr period is close to this va lue); (b) it is consisten t 
with the X-ray period derived bv lPrestwich et ail {2007), 
whereas the 35.7 hr period is 2a away from this value; 
and (c) it is surrounded by the lowest values of x 2 m 
Figure 0] (yet the value of x 2 begins to rise dramatically 
at periods just larger than 35.7 hr). The fact that we 
cannot completely discern between the two periods is the 
result of a half-cycle ambiguity between our Nov. and 
Dec. observations. Clouds precluded a second exposure 
on Dec. 12 which would have almost certainly eliminated 
this uncertainty. 
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TABLE 2 
Derived Black Hole Mass (Mq) 



Inclination 
(deg) 


Wolf-Rayet Mass 


(Ms) 


17 


25 


35 


90 


23.1 ±2.1 


27.7 ± 2.3 


32.7 ± 2.6 


78 


23.9 ±2.1 


28.6 ±2.4 


33.8 ±2.8 


65 a 


27.1 ±2.5 


32.3 ±2.8 


37.9 ±3.2 



a If the mass of [MAC92] 17A is 35 M , inclinations less 
than ~78° will not yield X-ray eclipses. 

4.2. Determination of the Black Hole Mass 

Given a radial-velocity curve of the secondary star, 
two parameters are necessary to calculate the minimum 
mass of the companion BH: the period (P) and the semi- 
amplitude {K^). The standard method of measuring stel- 
lar masses in binary systems is to calculate the m ass 
function, / (M) (e.g., iMcClintock fc Remillardl[2006l) . 



/(Mi) 



PKj 
2nG 



M x sin 3 % 



(1) 



where i is the inclination of the binary's orbit and q = 
M2/M1; here, Mi is the mass of the putative BH and Mi 
is the mass of the WR donor. Physically, / (Mi) is the 
smallest possible BH mass; it is equal to the BH mass 
only if i = 90° and q = 0. 

Using our measured values of P and K2 (see i j4.ll 
above), we calculated the mass function for the system. 
Finally, we assumed a range of reasonable masses for the 
WR star, which then allowed us to calculate a range of 
possible BH masses. 

Adopting our calculated period and semi-amplitude, 
the resulting mass function u sing Equation[T]i s 7.64± 1.26 
M Q , consistent with that of IPrestwich et all (|2007l ) (7.8 
Mq). If we assume a mass for the WR companion and 
an inclination for the orbit, then we can calculate a mass 
for the BH. 

The spectroscopic mass of the WR star was calculated 
by IClark fc Crowtherl (|2004l ) to be 35 Mq, but there is 
much uncertainty in this value. IPrestwich et al.1 ()2007l ) 
point out the small possibility that [MAC92] 17A could 
have a mass as low as 17 Mq. They also note that since 
this system is observed to have deep X-ray eclipses, the 
inclination of the system should be close to 90°. 

Table [2] shows possible BH masses assuming a reason- 
able range of values for both the mass of [MAC92] 17A 
and the orbital inclination (similar to IPrestwich et al.l 



120071 Table 1). Given these ranges and our calculated 
values for the period and semi-amplitude of the radial- 
velocity curve, the BH in IC 10 X-l must have a mass 
of at least 23.1 Mq, and it very well might be 32.7 M Q 
or larger. In any case, this is much l arger than the mas s 
of the BH in M33, 15.65 ± 1.45 M ijOrosz et al.ll2007l ). 
arguably the most massive previously identified stellar- 
mass BH, which was a nnounced very shortly before the 
BH mass of IC 10 X-l (jPrestwich et al.ll2007h . 

Note that if we had used the period of 35.7 hr men- 
tioned in and its associated semi-amplitude of 390 
km s , we would have derived a mass function of 
9.15 ± 1.45 Mq, yielding minimum BH masses of ~25.5- 
35.8 M (given the previously assumed range of WR 
masses, 17-35 Mq). 

5. CONCLUSION 

In this paper, we present new spectra of the WR 
stars [MAC92] 17A and RSMV 2 in the nearby star- 
burst galaxy IC 10. [MAC92] 17A has been shown to 
be in a bi nary system with the variable X-ray source 
IC 10 X-l (jPrestwich et al.ll2007f ). From our Keck spec- 
tra of [MAC92] 17A, we have constructed a compelling 
radial-velocity curve. The measured orbital period and 
semi-amplitude of [MAC92] 17A imply that the mass of 
the BH companion is at least 23.1 Mq, and more likely 
~32.7 M Q . Thus, we have shown that IC 10 X-l is indeed 
a WR/BH binary containing the most massive known 
stellar-mass BH. 
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